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I. INTRODUCTION 

Charge fluctuations in mesoscopic devices are increas- 
ingly important as the devices become smaller. At low 
temperatures, the statistics of these fluctuations is de- 
termined by quantum effects. Attention has focused on 
the full distribution of probabilities P„ for the transfer 
of n charges from one part of a system to another — the 
so-called Full Coimting Statisticsii In fcrmionic systems, 
the quantum nature of the system, and how it is driven by 
external stimuli, manifests itself, even for non-interacting 
fcrmions, in the current-current correlation function and 
the higher order correlation functions which are becom- 
ing increasingly accessible to experiment i^i^ii^ 

Most theoretical work has concentrated on the sim- 
plest possible device, namely a tunnel junction between 
two ID leads^i^iii^i^iiaiiiii^ It has been shown that, if a 
bias pulse, V{t), is apphed across a junction with fixed 
transparency at low temperature, the tunneling processes 
induced by the pulse arc combinations of two elementary 
types of event called uni-directional and bi-directionalii. 
It is also known that the statistics of the transfer of 
charge induced by such pulses depend strongly on the 
driving protocol. In contrast to a general shape of ac 
pulse which leads to an indefinite number of electronic ex- 
citations in the leads, an optimal ac signal, which is com- 
posed of overlapping Lorcntzian pulses, has been found to 
excite a strictly finite number of excitations per cycle and 
to bring the noise down to dc levelsi^ The minimal ex- 
citation states (MES)r^ created by such optimal pulses, 
offer the prospect of being able to generate signals with 
well-defined charge transfer down to the level of single 
electron emission. Coupled with the high Fermi velocity 
in electronic systems there is the prospect of rapid solid- 
state information transfer at a level useful for quantum 
information processingF^^iiiii^ii^iii 

Another area where the FCS have been studied is that 
of quantum pumps. These can lead to the pumping of 
electrons from one side of a tunneling barrier to the other. 
Several schemes for operating a tunnel junction, which 



can lead to the transfer of charge a^^i^^ and produce en- 
tangled electron-hole pairs in separate leads, have been 
proposed. Samuelsson and Biittiker have proposed an 
orbital-entangler, which works with quantum Hall edge 
statesii^ Accurate control of the transparency may al- 
low the generation and control of flying qubits, while 
Beenakker et al. have shown that such an electronic cn- 
tangler based on a biased point contact could reach the 
theoretically maximum efficiency of 50%<^ 

If the quantum effects are not to be obscured by ther- 
mal noise, the temperature 1 //? must be low enough that 
tf < (3, where t/ is the measurement time (or inverse 
repeat frequency for an ac measurement). Working at 
temperatures around lOmK would require operating at 
frequencies around 200MIIz, and this is the temperature 
and frequency regime used in some experimentsi^ How- 
ever, even at zero temperature the so-called equilibrium 
noise is present and diverges logarithmically with the in- 
verse repeat frequency or measurement time tj. This 
equilibrium noise is present in both the proposed MES 
protocol for generating charge transfer and the protocol 
for the optimal electronic entangler. 

Here we develop our approach to calculating the FCS 
for a tunnel junction^^ and examine protocols which we 
proposed for suppressing the equilibrium noise both in 
the case of the charge source and of the entangler^ We 
show how to solve for the FCS in the general case of 
fully time-dependent barrier profile with dynamic bias 
pulses applied between the leads. We compute the result- 
ing FCS and induced entanglement entropy for a number 
of profiles and; in particular, those close to optimal (in 
the sense that they have low noise in the case of elec- 
tron sources or maximum entanglement). Our approach 
is motivated partly by the result of Abanov and Ivanov 
( AI) who on quite general grounds deduced constraints 
on the analytic properties of the characteristic function 
(generating function for the probability distribution P„). 
Although they have recently argued that at any temper- 
ature, the counting statistics can be regarded as general- 
ized binomial statistics in which electrons scatter off the 
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barrier with some effective transparency independently,'^ 
our results are all for zero temperature. 

The paper is organized as follows. In Sec. |TT]we review 
FCS and we give a short alternative derivation of the 
AI formula. We then discuss the two 'standard' special 
cases: the biased contact with fixed transparency, and 
a contact with modulated barrier transparency. We es- 
tablish a mapping between these two cases and use this 
to simplify the derivation of the FCS for these two cases 
and to explore the relation to the Fermi Edge Singularity 
(FES) problem. In Sec. Illll we describe a general purpose 
numerical procedure to solve for the general case inacces- 
sible to analytical techniques. We use this to compute the 
effects of deviations from the ideal voltage pulses, which 
lead to minimal noise in the charge transferred across a 
tunnel barrier when operated as an electron source, and 
to study protocols close to optimal for the generation of 
electron entanglement. Concluding remarks can be found 
in Sec. M\ 



II. FULL COUNTING STATISTICS 

We consider a quantum point contact (QPC) at zero 
temperature with time-dependent transparency, T{t), 
connecting two single channel ballistic conductors, as il- 
lustrated in Fig. [1] We assume there is no inelastic 
scattering inside the QPC or the leads. The two leads 
(assumed identical) are disconnected initially and con- 
tain non-interacting electrons in their respective ground- 
states |0). Electrons in the disconnected leads are de- 
scribed by the Hamiltonian Hq = J2j ^j^l^j^ where we 
assume the system has a discrete energy spectrum. The 
electron creation operators have been written as a vec- 
tor c^^ = {c\^, c\i), where ^ is the creation operator for 
states in the left and right lead respectively. The ground- 
state of a single lead |0) is the Fermi sea, |0) ~ He <p ^] I)' 
where |) is the "true" vacuum and fi is the Fermi energy. 
Where necessary, we will assume a cut-off of order the 
Fermi energy. 




FIG. 1: (Color online) A quantum point contact with time- 
dependent bias voltage V{t) applied on the right lead. The 
transparency of the QPC, T{t), is controlled by the gate volt- 
age Vait)^ 

The two Fermi seas are initially uncoupled. Usually it 
is assumed^ that the tunneling barrier is lowered at time 
t = 0, allowing electrons to tunnel between the two leads. 



and is restored to fully reflective after a measurement 
time, tf. In general, the evolution of outgoing states 
should be described by solving the fully time-dependent 
Hamiltonian H{t) = Hq + H'{t), with 



i7'(t) = ^c]M(t,e„£,v)c, 



(1) 



The matrix M{t) = for f < and t > tf. However, if 
the scattering potential varies slowly on the scale of the 
Wigner delay time <C t\y ^ "^^^ffi with h = e ~ 1, 
the properties of the system can be determined from the 
instantaneous value of scattering matrix S{t) evaluated 
on states at the Fermi energy 



S{t,E = ii) 



B{t) A{t) 
-A*{t) B*{t) 



(2) 



Here A{t) and B{t) arc time-dependent transmission and 
reflection amplitudes and are determined by both the 
QPC gate voltage and the bias voltage applied between 
the leads. This relates eigenstates Cj of i/o, which we 
separate into incoming aj and outgoing states 6j, via 



Sit) 



(3) 



We arc interested in the distribution P„, which is the 
probability that there is a net transfer of n charges from 
the left to the right lead during the measurement period 
< t < tf. A convenient way of characterizing the Full 
Counting Statistics (FCS), P„, is via the function xW- 



X(A) 
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(4) 



The current, noise and higher order cumulants, ((Q'")), 
can be computed from x(A): ((Q")) = ^o(j0r\x=o, 
where m is the order of the cumulant. The formula for 
the FCS isjii^i^e 



X(A) = det(H-n(S'^e''^^S'e- 



1)), 



(5) 



where n is the number operator for the fermions. The 



matrix L projects onto states in the left lead: 







lead space. The matrix inside the determinant is infi- 
nite dimensional in the energy or time domain and 2x2 
in lead space. Because of the infinite dimensionality of 
the space of states in ([5]), careful regularization of the 
formula is required For example, at very high en- 
ergies £ — > -f oo, n ~ 0, the argument of the determinant 
approaches the identity and the contributions remain fi- 
nite and computable. However, when e —^ — oo, n = 1, 
and the matrix has the asymptotic form e^^^ Se~'^^^ , 
which makes the determinant ill-defined for an infinitely 
deep Fermi sea. 

A simple and correctly regularized approach to the 
computation of xW works with the density matrix for 
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the outgoing states in one of the leads. The density ma- 
trices for incoming states in both leads can be written 



Tij = 6{ii — Ej), which is the Fermi 



distribution function at zero temperature. Fourier trans- 
formed to the time domain, the density matrix has the 
form / d£6l(^ -' e)e'(*-*')(^~'') = ^ t-t'+tO ' 

The density matrix of the outgoing states in, say, the left 
lead n°'^*{t,t') can be obtained from ([3]): 

n™*(t,t') = {b[{t)bL{t')) 

= B*{t)nT{t,t')B{t')+A*{t)n'Sit,t')A{m 

In the second equation, we have used the fact that terms 
like {a[^afi) are zero, because the incoming states between 
different leads are uncorrelated. 

From it is possible to compute the cumulants 

directly. For example, the second cumulant (noise) is^ 

.2//.«,V.-,M')[l-..»(«'.«,] 

[\Amf\A{t')f + A-{t)B{t)A{t')B-lt')] (7) 

In the case, when the barrier transparency is switched 
on and off with A{t) ^ Aq = const for < t < tf and 
A{t) = otherwise, and with no bias voltage applied 
between the leads, the so-called equilibrium noise is ob- 
tained from the integral in ([7]): ((Q^)) = ^ \ogtf£_. Here 
^ is the commonly used ultraviolet energy cut-off of the 
order of Fermi energy. The logarithmic term is present 
for almost all profiles and not just abrupt switching. It 
was found, for example, for the case of a Gaussian switch- 
ing profile.— 

The eigenvalues, rij , of allow for the direct compu- 
tation of the FCS. We consider first a simple case, where 
only one eigenvalue changes to rij and all other eigen- 
values are unchanged. Since all eigenvalues initially are 
either or 1, we need only to consider a change — *■ nj 
or 1 ^ rij. In the first case, the probability of one ad- 
ditional particle being transferred into state j from the 
right lead is rij, while the probability, that no particle 
is added, (1 — rij). The counting statistics follow from 
(|4]) and are given by XjW = 1 ^ + e^^n-j, with the 
average charge transfer given by (Q) = 
If the occupation changes from 1 



In X I 



'J' 



a single hole 

is transferred with probability 1 — Uj while no charge 
transfer takes place with probability rij. For this case, 
XjW = rij 
fer (Q) = n 
be written 



(l — rij), with the average charge trans- 
1. The results for the two cases can both 



X,(A) = e*^(<'?>-"^)[l + (e*^-l)n,] 



(8) 



Since we work in the basis where is diagonal, the 
result for xW for the general case is simply a product 
over the factors, XjW- Taking account of all possible 
processes, we arrive at the formul a^^'^'^i^^'^^ , 



X(A) 



JMQ) 



N 

n- 



-iXrij 



[1 + (e'^ - l)n,] 



(9) 



It is correctly regularized as states unaffected by the per- 
turbation contribute a factor 1 to xW- It is also well 
suited to direct numerical calculation. 

In the following sections, we will discuss two special 
cases where the FCS can be obtained analytically. We 
rederive the known results for these cases by working di- 
rectly with the density matrix, Then we show how 
the FCS for the general case are easy to obtain by diago- 
nalizing n°"* numerically. To facilitate the interpretation 
of ([9|), we choose a specific form for the scattering ma- 
trix S{t), and assume that the transmission and reflec- 
tion amplitudes of the barrier A{t) and B{t) controlled 
by the QPC are real. If a bias voltage V{t) is apphed 
between the leads, its effect is to introduce an additional 
phase difference between the states in the two leads given 
by the Faraday flux ^{t) = | J* V(t')dt' . We incorporate 
this effect via a gauge transformation applied to states in 



the right lead: an 
matrix S{t) is: 

S{t)-- 



a^e 



#(t) 



The resulting scattering 



Bit) 



B{t) 



(10) 



A. Bias-voltage applied between the leads 

The case of a bias voltage pulse, V{t), applied across a 
barrier with flxed transmission amplitude, A, for < 
t < tf has been well studied.— i^i^ It has been shown 
that charge transfer processes at zero temperature are 
made up of combinations of two elementary events called 
uni-directional and bi-directionaL^^^ The uni-directional 
event relates to a single charge transfer process associated 
with the dc component of the bias voltage V{t). A single 
state is occupied on one side of the barrier and not on 
the other. This leads to the possibility of the transfer of 
charge across the barrier but only in one direction. Uni- 
directional events contribute to the average current as 
well as higher order cumulants. Bi-directional events are 
the consequence of the ac component of V{t), and relate 
to the excitation of equal numbers of particles and holes. 
These can both be transferred or reflected at the barrier 
so that charge can be transferred in either direction. No 
average current is generated in this case and only even 
cumulants are non-zero. The generic formula for the FCS 
for charge transfer isii 



e~^ + e-'^-2)]. 



X(A) = ll{R + Te^''')Yl[l + RT sin' ^( 
i=i j=i 

(11) 

T = |Ap is the barrier transparency and R = 1 — T. 
N.u and Nh correspond to the total number of uni- and 
bi-directional events, k = ±1 depending on the polar- 
ity of the voltage pulse. The angles aj/2 determine the 
probability of exciting a single particle-hole pair in a bi- 
directional event. N^, Nt and aj/2 can be computedii 
by diagonalizing matrix hh, where h and h are defined 
as /i = 2n - 1 and /i = UhU^i, with U{t) = e*'^(*). 
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We can understand the form of pT|) by considering 
the density matrix of outgoing states, which in the case 
of constant transparency has the form 

n°"*(t,t') = Rn'i^{t,t') +Te''^^*^n'^{t,t')e-''^''-*'l (12) 

We assume that the measurement time is short enough 
that we can ignore the equiUbrium noise contribution, 
which is a logarithmically divergent function of the mea- 
surement time The equilibrium noise is associated 
with fluctuations in the number of particles in the left or 
right lead and occurs even in the absence of an applied 
voltage. 

The dc component of the voltage pulse, associated with 
non-zero Faraday flux ip, generates additional occupied 
particle (or hole) states in the right lead when compared 
with the incoming states in the left lead. The corre- 
sponding particle (or hole), after impinging on the bar- 
rier, will tunnel across with probability T = yl^ or be 
reflected with probability R = 1 — T. This gives rise 
to so-called uni-directional events. The eigenvalue of the 
density matrix for outgoing states in the left lead is then 
rij = T, with an average charge transfer from right to 
left of (Q) = T \{ j relates to a state above the Fermi 
energy, or rij = 1 — T and {Q) = —Tiij relates to a 
state below the Fermi energy. Inserting this in ([9]), gives 
X„(A) = R + Te"^'^, where k = ±1 is determined by the 
type of transferred charge (particle or hole). 

An example of unidirectional events is provided by 
the so-called minimal excitation states (MES). These 
excite a number of particles (or holes) with minimum 
noisci^ii^ The corresponding voltage pulse V{t) is a sum 
of Lorentzian pulses: V{t) 

and \Tj\ determine the center and width of the jth pulse 
respectively. The unitary transformation in (jlOp is 



Ef (t-^+r- ' ^here t, 



N 



(13) 



Choosing the signs of the Tj to be the same leads to 
N unidirectional events. For N = 1, the pulse in (fT5|) 
generates a single uni-directional event with one parti- 
cle (electron or hole, depending on the polarity of the 
pulse) passing through the barrier with probability T 
giving x(A) = R + TeJ''^. For 2, with the polar- 

ities of the two pulses (set by the signs of the tj) the 
same, the two uni-directional events are combined and 
x(A) = (i? + Te^'^^Y irrespective of the relative widths 
(|rj|) or positions [tj). 

The ac component of the voltage pulse, associated 
with zero total Faraday flux, gives rise to the so-called 
bi-directional events — the simplest example of which is 
given by a pulse of the type (fT3|) with N = 2 and riT2 < 0. 
The operator hh characterizes the differences between the 
states in the two leads following the application of the 
pulse. All states, which are either both occupied or both 
empty in the two leads and which therefore contribute a 
factor 1 to x(A), are eigcn states of hh with eigenvalue 



1 {aj = 0). All other eigenvalues occur in pairs and are 
equal to e^*"^ F^i In the basis of the unperturbed states 
of either of the leads, the state described by h contains 
an admixture of the unperturbed state and independent 
particle and hole excitations in each of the 2 dimensional 
subspaces of the basis (labeled by j), in which h is block 
diagonal. Its eigenvalues are e**"^ . Each of these will 
lead to eigenvalues in n°"* of rij and (1 — rij), one as- 
sociated with the particle excitation and one with the 
hole. 

The contribution to x(A) from each of these bi- 
directional events (corresponding to the different j) is 
the product over the two factors of the type in ([8]), one 
with eigenvalue rij ((Q) — rij) and one with eigenvalue 
1 - n, HQ) = n, - 1): XbW = (l + (e'^ - l)rr,) x 
(1 + (e-*^ - l){l-nj)) ^ l + nj{l-nj){e^^ + e-*^-2). 
To make a connection between rij and rotation an- 



gle Q;j/2, we use the known resulljii h 

e-' 



1 

1 



and 







in the eigcnbasis of hh. Substituting 



h and h into p2p and diagonalizing n 
obtain 



± 



1 - ART sin^^ and 



nj(l - rij) = RT sin^ ^. 



explicitly, we 



(14) 



After taking the product over all events labeled by 
J, and adding in the contribution of the uni-directional 
events, we recover (|lip . For a given voltage pulse between 
the leads and corresponding unitary transformation U (t) , 
the aj/2 can be thought of as the rotation angles of the 
ground state associated with U (t) and are found by di- 
agonalizing the hhM^ The rotated state is an admix- 
ture of the original state, with probability cos^ and 
the state with an added particle and hole, with probabil- 



ity sin^ ^. The factor 1 + RT sin'' ^{e' 



2) 



is the weighted average of the result for the unper- 
turbed state (contribution 1 with weight cos^ ^) and for 
the state with an added particle and hole (contribution 
{R + Te'^){R + Te-'^) with weight sin^ ^)^ 



B. Barrier with modulated Transparency 

Another case for which results in closed form have 
been reported is that of a time-dependent barrier be- 
tween two leads at the same chemical potentia L^^i^^ 
The problem can be mapped onto a special case of a 
voltage biased time-independent barrier with constant 
transmission and reflection amplitudes This mapping 
becomes explicit once the problem is approached via 
the density matrix of the outgoing states. In the ab- 
sence of a bias the scattering matrix in simplifies 

I and the density matrix of 
-A{t) B{t) J 

outgoing states becomes n°''\t,t') = B{t)n{t,t')B{t') + 
A{t)n{t,t')A{t'). Introducing e^-^^*) = B{t) + iA{t) (wc 



to S{t) 



are still assuming that A and B are real), we insert 
g^Ht) into and eliminate A and B. We obtain 

n™* = i(e*'^n"e-''^ + e-^^n^e^"^). Here we have used 
the fact that = nj^* = n™. 

Since a unitary transformation on does not affect 
its eigenvalues, nj, we can study 

e-"^n°"*(f,t')e"^ = in*"(t,t') + ^e2"^(*)n™(i,i')e"^*'^^*'^ 

(15) 

The relations and have the same structure. The 
FCS of a system with modulated barrier transparency 
without a bias between the leads are therefore equivalent 
to those for a system with bias voltage applied across 
a barrier with constant transmission and reflection am- 
plitudes A = B = and Faraday flux ip — 20. The 
FCS for a modulated barrier transparency can therefore 
be obtained from ([TT|) . In addition, other concepts de- 
veloped to understand the bias voltage case carry over to 
the modulated barrier case. These include the geomet- 
rical interpretation of the FCS^^ as well as the MES^^ 
which, when implemented as a modulation profile of the 
barrier leads to the optimal entanglcr of electron hole 
pairs 

To calculate xW for the case of the modulated barrier 
from (|lip . we need, as before, to diagonalize the matrix 
hh = /ie^"^/ie~^*'^ and compute the angles aj/2 from its 
eigenvalues. Since no bias voltage is applied between the 
leads, the system is completely symmetric in lead space. 
As a result only bi-directional events (leading to no net 
average charge transfer) can occur and all the eigenvalues 
of /ie^*'^ft.e~^"^ come in pairs. The characteristic function 
x(A) is then given by 

= n(l + ^Bin2^(e^^ + e-^-2)) 

= n(l-Bin^f sin^^), (16) 

where N is the number of paired eigenvalues determined 
from he^"^he-'^'*. 

To illustrate the value of this mapping, consider the 
case in which the transparency of a barrier is subjected 
to a sinusoidal modulation: A{t) ~ sinujt and B{t) ~ 
cosojt. If the total number of cycles is large enough, the 
contribution from the equilibrium noise (proportional to 
logtf^) can be neglected. This problem can be mapped 
to a case of a barrier with constant bias voltage with 
transparency T = 1/2 and -0 = /* V(t')dt' = 20 = 2ujt 
corresponding to a constant dc bias with V ~ 2u!. 

To obtain the FCS, we need to diagonalize matrix 
^g2ia;t^g-2iut ^ i-j^jg corresponds to a constant dc 
bias problem involving two uni-directional events per pe- 
riod, (the phase changes by Att per period), there are 
two eigenvalues different from 1 and both are equal to 
-1. The polarity (particle or hole) of transferred charge 
can be inferred from the requirement that there is no net 



FIG. 2: (Color online) Four possible outcomes per period 
from a sinusoidally modulated barrier transparency. AH four 
outcomes are equally likely. The particle-hole pair can be on 
either side of the barrier or as an entangled pair of particle 
on one side and hole on the other. 



average charge transfer. Hence wc conclude that the si- 
nusoidally modulated transparency case is equivalent to 
one electron and one hole impinging in a single period on 
a barrier with transmission T = 1/2. The two processes 
are independent since both correspond to uni-directional 
events in the equivalent bias voltage problem. The FCS 
for constant bias voltage case are known to be given by 

Xo{V,X) = {R + Te"^^)-i^, with k = ±1 depending on 
the polarity of the applied voltagc^^*^ The corresponding 
FCS for the sinusoidally modulated transparency case is 
a combination of two factors xo one with k = 1 and one 
with K = —1. This gives 

x(A) = xo(^, A)xo(w,-A) 

= {^f, (17) 

which is a result previously obtained by Andreev and 
Kamcncv using the Keldysh formalismjAS A barrier oper- 
ated in this way acts as a quantum pump which excites 
exactly one electron-hole pair per period, provided the 
logarithmically divergent equilibrium noise is neglected. 
The four possible outcomes per cycle arc shown in Fig. 

121 

The result p?|) is the ac version of the barrier profile re- 
quired to generate optimal electron entanglement.— The 
profile 

^ B{t) + iA{t) = ^—^ (18) 

t + IT 

generates the FCS x(A) = (1 + cosA)/2. The use of 
a quantized Lorentzian pulse of the type ^TE\i has the 
advantage that the equilibrium noise is strictly absent 
in (fT5|) . Although the entanglement generated by (fTSl) 
is between electrons and holes and therefore not useful 
in quantum computation due to the charge conservation 
laws^S, the protocol would be useful if combined with 
another degree of freedom (either spin or orbital) because 
it operates in a one-shot mode. The difficulty with its 
operation is associated with the precise generation of the 
actual profile. We explore the effect of possible errors in 
its experimental implementation in Scc lIII Bl 
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C. Fermi Edge Singularity 

We discuss here the relation between the FCS of a con- 
tact subjected to sharp bias voltage pulses and the Fermi 
Edge Singularity (FES). We map to the problem of an un- 
biased contact for which the result for the FCS is known. 

We consider two delta pulses with opposite signs sep- 
arated hy tf applied to one lead. The contact has fixed 
transparency T and reflectance R. The corresponding 
voltage profile is 

Vit) = ^m-6it-tf)], (19) 

with ip = const. This pulse induces a phase shift tp on 
the incoming states of, say, the left lead (measured with 
respect to the right one) within time window < t < tf 
and zero phase shift otherwise. For t/i 3> 27r and large 
tf^, the noise can be written as an expansion in i£.tf)~^ 

((02))=2i?T[-|sin2|lnVe+|^] + -- - . (20) 

The derivation of (pO|) is essentially the same as that given 
by Lee and Levitov (LL) for the quantum current fluc- 
tuations induced by a magnetic field in a metallic loop 
containing a QPCi^ The FCS for the problem we are 
considering have been computed^ but only for the case 
that the transparency of the barrier is low. In this section 
we address the general problem with arbitrary contact 
transparency. 

As the two electrodes are equivalent, only bi- 
directional events occur in this problem. By (jlip . the 
counting statistics are 

x(A)=n[l-sin^f Bin^^], (21) 

j 

where we have defined sin -I = 2 V RT sin ^ . The cen- 
tral problem is how to compute angles Eq. (PT|) 
has the same form as ^TE\\ . which is that for a bar- 
rier with time-dependent reflection and transmission am- 
plitudes B{t) + iA{t) = e*'^(*)/^ For < t < tf, 
A{t) = sin^ ^ const and A(t) = if t < or t > t/. One 
interesting complication in this problem is that the FCS 
are only well-defined if the limit of increasingly sharp 
voltage profiles, V{t), leading to the delta functions in 
(fT9)) . is specified. This is because the shape of the volt- 
age pulse, leading to the Faraday flux change from zero 
to ip, affects the number of excitations introduced by the 
switching process. In the mapped barrier opening prob- 
lem, the precise time-dependence of the (rapid) opening 
and closing of the barrier matters. An example of a pos- 
sible profile for this opening of the barrier is shown Fig. 

El 

For the case that the barrier switching time t is very 
short compared to the measurement time tf, we argue 
that the switching processes, which excites high energy 
(~ 1/t) particle-hole pairs, should not interfere with the 
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FIG. 3: (Color online) Illustrative transmission amplitude 
A{t) as a function of time t for the the pulsed lead case af- 
ter mapping to the equivalent barrier opening problem. The 
FCS are only well-defined if the exact time-dependence, whose 
limit is the delta function in (|19|l , is specified. Here we show a 
case where the amplitude, A{t), actually overshoots the value 
sin(-(/)/2). We argue that the FCS should be well approxi- 
mated by a product over two independent contributions from 
the two regions denoted by 1 and 2 (shaded region). 

long time measurement process, which leads to the low 
energy excitations, predominantly on the energy scale 
~ 1/tf, expected for a Fermi Edge Singularity problem. 
We make the Ansatz that the FCS can be written as a 
product of the types of processes as 

x(A) = Xi(A)x2(A), (22) 

with X2 giving the FES contributions (which are associ- 
ated with the shaded region, shown schematically as the 
shaded area in Fig. [3]). xi gives the contribution associ- 
ated with the opening profile at t = and tf. We focus 
here on the calculation of X2 and postpone the computa- 
tion of xi to Sec. nil Bl where numerical techniques are 
adopted. 

The FCS of an unbiased barrier, for which the reflec- 
tion amplitude is abruptly changed from zero to the con- 
stant value sin ^ for a duration tf, is known and was com- 
puted using the bosonization technique.— The result can 
also be found by solving a Ricmann-Hilbert problem^ 
valid even at non-zero temperature in Appendix IB] The 
FCS is 

X2(A)=exp(-A2G), (23) 
with G and A.^ being given by 

A* V A , tP \ , , 

sin — — sin — sin — = 2v i?T sin — sin — (24) 
2 2 2 2 2 ^^ 

G = ^Ini/C- (25) 

The logarithmic terms in are connected with the 
Fermi Edge Singularity found in metals.— It is in- 
teresting to note that the form ([M)) includes the well- 
known result for the Anderson orthogonality catastro- 
phe problem in a single lead^^ as a special case. The 
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quantity of interest is the overlap (0|0') between |0), the 
ground state of an unperturbed metal and |0'), the wave- 
function of the same metal at a time tf after switch- 
ing on a (core hole) potential. In the simplest case, 
the potential can be well described by a single energy- 
independent phase shift which is exactly equivalent to 
the phase ip, with e^"^^^ = B + iA associated with the 
modulated barrier we have been considering. The over- 
lap can be written as (0|0') = (Ole''^^*) |0) = HjCos^, 
where the ^ are the eigenvalues e*"^ of he^^^'^he"^'^/'^ . 

This is equal to ^x2(A = 7r,r= \). With A = tt, T = i, 
A* = V' and we recover Anderson's result^: (0|e*'''^*) |0) = 

III. GENERAL CASE 

In the general case, which is equivalent to having both 
a bias voltage between the leads and a time-dependent 
profile for the barrier, we are not aware of the existence 
of a simple relation mapping the problem onto an equiv- 
alent bias voltage problem. However, ^ allows for the 
calculation of the FCS as long as the spectrum of the 
density matrix of outgoing states is available.— We can 
diagonalizc n°"* in ([6]) numerically to find the eigenvalues 
Uj and use © to compute the FCS. 

Numerically, it is convenient to work in the energy do- 
main with a discrete energy spectrum. (An alternative 
is to compute the dynamics of the system directly. For 
a finite system the spectrum is then automatically trun- 
cated and the determinant is properly rcgularize d'^^'^^ .) 
We introduce a periodic boundary condition in time, with 
period tp, discretizing the energy spectrum of the scat- 
tering matrix with an energy separation = 2-k /tp. By 
choosing tp tf we can compute the counting statis- 
tics with large number of cycles, giving the characteristic 
fimction as x(A) w [xo(A)]*^/*^ where xo(A) is the FCS 
for single period. We can also set tp ^ tf and study the 
behavior of a device when operated in one-shot mode. 

Fourier transformed into the energy domain, the in- 
dividual matrix elements of A{t) and B{t) are Xm„ ~ 
j_ J-Jpe-'(™-")'^o*A(t)dtwhereAstandsforAorB. The 
neighboring rows of X„in have the same elements though 
they are shifted from each other by one column. If X{t) is 
sufficiently smooth, by which we mean its Fourier trans- 
form decays faster than ojq, with v < —\. We can cut 
off the Fourier series and limit the approximated sum- 
mation within |m — n| < M with M chosen large enough 
to achieve the desired accuracy. The truncated matrices 
A and B are blockwise tridiagonal with each block size 
M X M . After some manipulation, the diagonalization 
of the infinite dimensional matrix is approximately 
equivalent to the diagonalization of a 2M dimensional 
matrix in energy space. Details of this procedure are 
summarized in App. [Xj 

In this section, we look at cases where the direct di- 



agonalization of the matrix n°"* allows the study of the 
FCS of a tunnel barrier with time-dependent scattering 
amplitudes with a bias voltage applied between the leads. 
We concentrate, in particular, on the example of a bar- 
rier used as a low noise source and study how the noise 
and degree of entanglement are affected by deviations 
from the optimal pulses which have been proposcdj ^^i'^^ 
Mostly, we will study barrier and voltage profiles which 
are cither combinations (or close to combinations) of the 
quantized Lorentzian pulses. These can be applied as 
voltage pulses to a lead with describing the corre- 
sponding Faraday flux e*'^'^*-', or with describing the 
barrier profile e*'^^*) = Bit) + iA{t). 

A. Quantized Pulses 

To operate a tunneling barrier as a single electron 
source at low temperatures requires a voltage pulse which 
excites a single electron excitation in one lead. This can 
be achieved by creating a MES with a single Lorentzian 
pulse applied between the leads. At low temperatures, 
the noise produced by such a device comes from two 
sources: shot and equilibrium noise The shot noise for 
the simplest case of a barrier with constant transmission 
probability, T, is proportional to T{\ — T) for a single 
MES pulse. This suggests that one should aim to open 
the barrier fully {T — 1) to increase the chance of sin- 
gle electron emission. However, if we open the barrier in 
an arbitrary way, the equilibrium noise becomes impor- 
tant. A solution^^ is to combine the creation of the MES, 
choosing a bias voltage pulse giving Faraday flux 

e''>(*) = (26) 
t - <o + ^T^i 

in the incoming states, with a carefully-chosen opening 
profile for the barrier which minimizes the total noise. 

Since the profile generating the maximal entanglement 
(fT8)) . as illustrated in Fig. IDJa), fully opens the bar- 
rier twice without any equilibrium noise contribution de- 
pending logarithmically on tf, a first guess is that it 
might also work as a possible profile for the opening of 
the barrier when used as a single electron source. How- 
ever, the barrier opening scheme based on (jlSp is not 
ideal for a single electron emission. Though the log- 
arithmic term in the equilibrium noise is absent, the 
profile docs generate background noise with (sec I16p 

((Q'»o = = Ef ^in'f < f . For the pro- 

file mi) and iV = 1 and a = tt we find ((Q^)) ^ i which 
is actually the maximum for a single pulse. 
We consider instead a pair of such pulses: 

g»0(t) ^ t-tl/2-iTlt + tl/2 + iT2 

t-tl/2 + iTlt + tl/2-iT2' ' 

where the separation of the two pulses (with widths ri 
and T2 respectively) is ti. A little algebra shows that 
A{t) and B{t) do not change sign only if ri = T2 and 
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FIG. 4: (Color online) Transmission coefficients A{t) for some 
opening profiles for a barrier. The centers of the pulses are 
at t = 0, with time scaled in units of the pulse width r. (a) 
Single Lorentzian pulse (see 118(1 for optimal entangler. (b) 
Two Lorentzian pulses, with opposite polarities (see [27} with 
Ti — T2 — T and pulse separation ti/r ~ 0.6. (c) Same as in 
(b) with separation ti/r ~ 0.83 (optimal), (d) Same as in (b) 
with ti/r ~ 1.2. The inset shows the corresponding reflection 
amplitude _B as a function of t/r. Only in cases (b) and (c) 
do both A and B remain positive during the pulse lifetime. 



Ti/ti > 5 + • Corresponding typical profiles are shown 
in Fig. SUb-c). The ratio n/ti = 5 + gives a trans- 
parency of the barrier which has a single maximum with 
A = 1 at i = 0, as shown in Fig. IH^c). We have found em- 
pirically that this separation gives the best combination 
(low total noise and highest probability for the transfer of 
one electron). Exciting a single excitation in the incom- 
ing states of one lead and coordinating the timing of this 
excitation with the opening of the barrier, allows the sin- 
gle particle excitation in the incoming states to impinge 
upon the barrier when it is fully open. An advantage 
of the profile ([77]) is that, owing to cancelation between 
the two components (at t = — <i/2 and t = ti/2) at long 
times, the transmission amplitude A{t) approaches zero 
faster than for the profile p^ . In addition, the noise 
generated by this opening profile ((Q^))o ~ 0.23, which 
is less than half that generated by . We expect that 
the profile ((27| . with an optimally chosen ratio for ti/ri, 
should be a good candidate for designing an on-demand 
single electron source at ultralow temperatures. 

For Ti/ti = I + ^Tf ^'^'^ with a single pulse e*''' = 
t-to-iTo applied to the lead, wc have diagonalizcd the 
density matrix nout numerically. Fig. EJa) shows results 
for the noise ((Q^)) in the system as a function of the 
separation to between the center of the MES pulse and 
the maximum of barrier transparency (for which A ~ 1 
at t = 0). The maximum values of the noise occurs when 
to ~ Ti and the transparency coefficient is almost 1/2. 
This is the regime where the barrier is acting as a 50% 
beam-splitter. The minimum of ((Q^)) corresponds to 
to = 0, when = 1. At to = the transferred charge 
as well as the quality factor {Q) / {{Q'^)) attains its max- 




FIG. 5: (Color online) Upper panel: Single unidirectional 
event in a quantum contact with tunable transparency (nu- 
merical calculations), (a) Transparency of the barrier as a 
function of time given by (|27|l . We also show an MES applied 
to a lead (|26p . (b) Average charge transfer (Q) between the 
two leads as a function of to. The maximum of (Q) is around 
0.96 for the narrowest MES considered ri/ro = 16. (c) Noise 
in a biased contact with transparency (|27p as a function of 
to/T\. The minimal value of the noise corresponds to to = 
and, for the narrowest MES pulse {ti/tq — 16), it almost 
drops to the value due to barrier modulation: {{Q^))o (dot- 
ted line) . The two maxima in the noise occur when to /n ~ 1 
and the transparency is almost 1/2. Lower panel: Probabil- 
ity distribution P„ for n/ro = 1,4 and 16 at to = 0, as well 
as the probability distribution generated by profile (|27fl only 
without the MES applied on the lead. We see that in the 
limit Ti/ro 00, P„ approaches Pn±i- (See text) 



imum. Fig. EJ^b) shows the average transferred charge 
(Q) as a function of to. If the pulse applied to the lead 
is narrow compared to the opening time of the barrier 
(tq ^ Ti), the minimum value of the noise is essentially 
set by the noise associated with the opening of the bar- 
rier, namely ((Q^))o- The transferred charge approaches 
1 when tq/ti —f and the probability distribution for 
transferred charge, P„, approaches Pn±i (the sign de- 
pends on the polarity of the voltage pulse applied to the 
incoming states in one of the leads) is shown in the lower 
panel in Fig (5] Here Pn±i is the probabifity distribution 
for charge transfer associated with the opening of the 
barrier without a bias pulse applied between the leads. 

Fig [5] shows that, when ti/tq = 1, the emission of 2 
electrons is strongly suppressed. If it was important to 
have a source which only emitted single electrons and 
there was a method for discarding 'non-events' in which 
no electron was emitted, this protocol would work well. 
The suppression of double electron emission can be un- 
derstood as follows. The width of pulse determines the 
energy profile of the excited particles or holes. When 
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FIG. 6: (color online) Evolution of the probability for 
electron-hole pair emission as a function of pulse separation 
ti/n for a barrier profile (|27|l with T2 = ti. The dashed 
line gives the probability that nothing happens (-Po)- The 
solid line is the probability of generating a single entangled 
electron- hole pair (Pi). The dotted line show the probabil- 
ity for two particle-hole pair emission (P2). The vertical cut 
is the position where the single electron source profile (|27|) 
operates. 



To ^ Ti, the particle excitation induced by the changing 
barrier profile has the same energy profile as that of the 
incoming excitation induced by the bias voltage applied 
between the leads. As the two pulses are coincident in 
time the Pauli exclusion principle leads to destructive in- 
terference between the two excitations in one lead and 
either only one particle or no particle will be transmitted 
through the QPC as a result. 

The profile (|27p does not involve a change of sign of 
either A{t) or B{t) and, consequently, should be easier 
to implement experimentally than the profile (fT5|) . We 
therefore consider how well the profile (j27p might work 
as a (non-optimal) electronic entangler. Fig. [5] shows 
the probability of electron-hole pair production Pn ver- 
sus pulse separation ti /ti . The vertical line corresponds 
toti/ri = 1/(5 + "^) ~ 0.83, beyond which a sign change 
in A(t) is necessary and the experimental implementation 
is expected to be more involved. We see that the emis- 
sion of two electron-hole pairs is very unlikely. The only 
two significant outcomes are: no excitation, created with 
probability Pq, and the creation of a single (entangled) 
electron hole pair which is created with probability Pi. 
The corresponding entanglement entropy at this point is 
S = PiSbcU ~ 0.23, ie just under half the theoretical 
maximum. 



B. Non-Quantized Pulses 

So far, we have discussed (combinations of) quantized 
Lorentzian pulses applied to a lead (see I13p and / or bar- 
rier opening profiles (see l27p that excite a bounded num- 
ber of electron-hole pairs without the accompanying equi- 
librium noise which grows logarithmically with to. Here, 



we focus on the increased noise resulting from deviations 
from the ideal quantized pulses. 

We consider the case of a modulated barrier, with the 
following non-quantized opening profile: 



B{t) + iA{t) = e*"^'*) 



t — ITl 

t + in 



71 



t - to - iT2 
t - to + iT2 



72 



(28) 

where 71,2 <= K-. We look at the simplest case 72 = 1 — 71 
and 72 = —71 and choose ti = T2 = t. In both cases 
to plays the role of the measurement time t/, provided 
to » T. The profile (|^ with to ^ t describes a barrier 
with a transmission amplitude. A, which changes from 
to sin 27r7 within a period of r around t = and closes 
at t = to. There will then be a contribution to the noise 
which increases logarithmically with to (this is just the 
equilibrium noise contribution). The mapping between 
the case of a barrier profile and a voltage bias across a 
junction with constant transparency means that this pro- 
file also models non-quantized voltage pulses applied be- 
tween the leads. This mapping includes a doubling of the 
total phase change 2(/) = -0 (see Sec. Ill Bp . which means 
that (pS)) describes two voltage pulses with Faraday flux 
47771^2. If the 7i are not both integer or half-integer, there 
is a net phase shift of 27771 and 27r(7i + 72) between the 
rotated and the unperturbed states with consequent FES 
effects. 

The FCS for the profile HH), which gives the opti- 
mum level of entanglement, are equivalent to those for 
two quantized Lorentzian pulses, when the problem is 
mapped to the case of a junction with a bias, because of 
the associated doubling of the phase 2(j> = ip, This sug- 
gests that the pulse (fT5)) applied to the barrier may be a 
special case of two separate pulses, each of which corre- 
sponds to a single quantized pulse in the case of a bias 
between the leads. In particular, a pulse 



i0(t) 



t — ITl 
t + ITl 



t-tn 



IT2 



t — to + iT2 



(29) 



consists of two pulses centered at and to with widths ri 
and T2. The absence of a logarithmic contribution to the 
noise in this case should be expected, because the barrier 
is closed at all times except for a period r around and to. 
The effect of pulses (|29|) on the states of the system can 
be found by working in the basis in which the scattering 
matrix in pH)) is diagonal. With A and B both real 
{ip — as there is no applied bias), the scattering matrix 
is diagonal for all times in the basis ci^2 = (cl ± icfi)/\/2 
and its eigenvalues as a function of time are given by 
the values of e**"^'*^ in ((29)) . The two components of the 
pulse centered on t = and t = to induce phase shifts of 
±7r in each channel and are individually as far as possible 
from quantized pulses. However their effect on the FCS 
of the charge transfer between the left and right leads 
depends only on the difference in phase shift between the 
two channels and each of the two components therefore 
contribute to the FCS as quantized pulses. 



10 




FIG. 7: (Color online) The noise produced in the charge trans- 
fer across an unbiased barrier when its profile is varied accord- 
ing to (|28[l . Results are shown as a function of the separa- 
tion, to/T, between the two non-quantized pulses aX t = and 
t — to, and of the exponent 71 with (a) 72 = 1 — 71 and (b) 
72 = —71- In case (a) the values 71 = 0.5 and 71 = 1 corre- 
spond to quantized pulses (when mapped to the problem of 
a biased lead, there is a doubling of the total phase change). 
The noise for all other values of 71 increases logarithmically 
with to. This is the equilibrium noise contribution. In case 
(b) the noise 71 = 0.5 and 71 = 1 saturates when the pulse 
separation to is much larger than the pulse widths r. For 
all other 71 there is again the equilibrium noise contribution 
which grows logarithmically with measurement time to. 



The noise generated by the opening profile ([28]) with 
Ti = T2 = T is shown in Fig. El^a) for the case 72 = 1 — 71 
and in Fig. [7{b) for the case 72 = —71. In both these 
cases the penalty for missing quantization of the pulses 
is small. The increase in noise at fixed to, as 71 deviates 
from integer or half-integer values, is slow. (For example, 
with 71 = 1/2 and to/r = 50, a 10% deviation in 71 in- 
troduces additional noise of only of 2% of the quantized 
value.) This suggests that any reasonable experimen- 
tal implementation of such pulses should allow the ex- 
ploitation we have described (as entanglers or as electron 
sources). Case (a) corresponds to two Lorentzian (non- 
quantized) pulses with the same polarity separated by a 
time tg. There are two peaks, the height of which grow 
logarithmically with to, with a flat valley in between. 
When 71 = 1/2 the profile is that of ([M)) . and we find, as 



expected, results equivalent to the profile ifTS]) , ie a total 
noise equal to 1/2 and independent of tg. For other val- 
ues of 71 the noise grows with a logarithmic dependence 
on to- When r < toi the problem is that of two equivalent 
FES problems: The effect of the first pulse is to give rise 
to scattering phase shifts of 27r7i and — 27r7i in the two 
independent channels in which the scattering matrix S{t) 
is diagonal. In case (b) there arc two oppositely polar- 
ized (non-quantized) Lorentzian pulses. At small pulse 
separation, to, the two pulses partially cancel and the 
noise is low. In this limit the barrier transparency, T(t) 
remains close to zero and vanishes when to = 0- Satura- 
tion of the noise at some finite value for large to/r occurs 
only when 71 = 1/2 and 1, when the two components of 
the pulse contribute to the FCS and each corresponds to 
quantized MES in the equivalent lead problem. For all 
other values of 71, we again find noise which grows with 
the logarithmic dependence with to expected for a FES. 

Another issue with the quantized Lorentzian pulses, 
used either as voltage pulses or to open the barrier, is 
that they are defined over an infinite time interval. In 
practice, the long tail behavior will be restricted to some 
closed interval [— <o,io] and this cut-off is likely to in- 
troduce additional noise. We have found, however, that 
amending the profile at points ±io and appending an 
exponential tail Ae{t) = A(to)e~'*~*"'/'^ to model this 
effect that there is virtually no difference between this 
trimmed profile and the ideal profile. This suggests that 
such deviations from ideal pulses arc unlikely to affect 
the operation of devices in this regime. 

Finally, we return to the problem introduced in Sec. 
Ill C[ where two delta pulses (sce [T9|) with opposite polar- 
ities were applied to one incoming channel. We argued 
that the FCS could be written as the product of two fac- 
tors (22). One factor, xi; is associated with the sudden 
jump in scattering phase shift and the other, X2, with 
the FES problem. To compute Xi(A), we need to specify 
the exact shape of the Dirac delta function. Here we use 
a Lorentzian in the limit of vanishing width r 04S 

The FCS are fully determined by the spectrum of n""*, 
or equivalently, the eigenvalues e*'"j of hh in the form 
sin^ ^ . We note that this problem is equivalent to a lim- 
iting case of the non-quantized barrier profile (|28p shown 
in Fig. El^b): After mapping, the barrier profile in (1^5]) 
has 71 = —72 = and r ^ 0. 

In the following, we show that the factorization of the 
counting statistics which we propose (|22[) is apparent in 
the separation of the values of sin^ ^ . We rewrite ^ = 
27rfc -I- 4>o, where fc S Z and -ipo G [0, 27r). When -ipo = 0, 
from pT|) the counting statistics are described by k bi- 
directional events with eigenvalues such that sin^ ^ = 1 
with j = l,...,k. Hence xi(A) = [1 - 4:RT sin^ ^]''. On 
the other hand, for states which contribute to (|23p. they 
are rotated with eigenvalues e^'"^ such that sin^ ^ is 
small (but total number of such eigenvalues is large) . For 
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FIG. 8: (Color online) (a) Evolution of the extra eigenvalue 
sin'^ ^ with ip in (|19p . (b) Noise computed from three ap- 
proaches: Numerics (dash dotted), ansatz (|32|) (solid) and 
LL's result (|20|l (dashed), (c) {{Q'^))x (see text) computed 
for LL's result and ansatz (|32p . 



the case when ipo 0, we have verified numerically that 
the k pairs of eigenvalues of hh with e^'"^ remain. In 
addition there is one pair of eigenvalues associated with 
tpo which gives sin^ ^ in the range [0, 1). Here ao/2 is 
the rotation angle associated with the phase i/jq. As a 
result, the total FCS ((2T|) should be well approximated 
by the form: 

X(A) = f](^l-4i?Tsin2^sin2^) 

« (l-4i?Tsin2^sin2-) 

• (^l-4i?rsin2 ■X2{X). (31) 

This value of ao can be computed numerically once for 
a simple case like k = 0, and tabulated over the range 
V-o G [0,27r). 

To test above idea, as well as to obtain the unknown 
eigenvalue sin^ we compute the FCS and the values 
sin^ ^ for (|30p numerically for different total phase shift 
-0 in (fT9|) . The pulse separation is fixed at to/r = 500 
and the contact transparency is T = 1/2. With cut-off 
energy ^ as the only fitting parameter, the ansatz ([3T|) 
agrees with the numerical result for xW i^i the parame- 
ter space {"0, A} everywhere with a maximum difference 
of 3%. This agreement justifies our decomposition of 
x(A) based on the separation of eigenvalues of The 
evolution of sin^ ^ with ijj is drawn in Fig. [5Ja) and 
echoes that reported for a tunnel junction driven by a si- 
nusoidal voltage as a function of the voltage amplitude4i 
As ip increases, the value sin^ ^ increases gradually and 
saturates at 1. 

The noise can be computed using (j3ip and is 

((Q2)) =2i?T[4sin2^1nio? + A: + sin2^]. (32) 



In Fig. [8jb) we show the noise computed from three ap- 
proaches: i) exactly computed numerically from ([9]); ii) 
using out ansatz (|32p and iii) from the expansion used by 
LL ([^D|) . The cut-off energy ^ is used as fitting parame- 
ter to give the best agreement with the exact result. The 
similarity between the noise produced here and the noise 
generation in Fig. Wijo) for fixed tj can be understood us- 
ing the mapping between the problem we are considering 
here of an applied bias voltage and that of a modulate 
barrier profile discussed in Sec. IIIBI In both cases the 
contribution proportional to the log tg^ is suppressed at 
quantized phase shift "0. 

The differences between the two approximate treat- 
ments and the exact results are associated with the 
"large" eigenvalues. We define the quantity {{Q^))x = 
{{Q^))nu — 2RTx{ip), with x{iIj) = ^ for the expansion 
(HOI) and x{i') ^ fc + sin^ ^ for our ansatz dST])- {{Q^))nu 
is the exact value of the noise computed numerically. Ac- 
cording to (PD|) as well as ([5^ . {{Q^))x is proportional to 
sin^ ^, and, if the cut-off energy ^ is independent of ip, 
is periodic. It should also vanish at = 0, 27r, 47r . . .. We 
show {{Q^))x for the two cases respectively in Fig. IDJc). 
We see that the ansatz (f32|) works qualitatively correctly 
while the expansion (j20p gives the wrong positions for the 
minima. However, we find that ((Q^))mi is not strictly 
periodic (the amplitude of the oscillation decreases). We 
attribute this to the fact that the separation into the two 
factors xi (corresponding to k large values of sin^ aj ) and 
X2 (only small values) is not complete and that the value 
sin^ afe_i_i/2 can be significant particularly for small k. 

IV. CONCLUSION 

We have discussed the Full Counting Statistics (FCS) 
of the charge transferred across a quantum point con- 
tact. We have illustrated the power of a mapping between 
the case of a biased barrier with constant transmission 
and reflection amplitudes and the case of a barrier with 
time-dependent profile but no bias. With this mapping 
we have showed that known results for the two cases, 
which had been previously obtained using different, and 
generally involved calculations, can be understood us- 
ing the basis of quantized Lorentzian pulses^ or Min- 
imal Excitation Statesi^ (MES). Examples include the 
optimal protocol for electron entanglement the FCS 
of a sinusoidally driven barrier^i^ and the Fermi Edge 
Singularityi^ 

For the purposes both of conceptual understanding and 
computation, we have argued that the problem is sim- 
plest when approached through the eigenvalues of the 
density matrix of the outgoing states in one of the leads. 
For the general case, which corresponds to applying both 
a bias and varying the barrier profile with time, we have 
developed a numerical scheme for computing exactly the 
FCS for a device and used this to compute the FCS for 
a tunnel barrier operated as an electron source. We have 
also studied how the deviation from an ideal pulse affect 
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the quality of operation of a device with low noise or an 
entangler. We showed that the noise levels were remark- 
ably insensitive to deviations from quantized Lorentzian 
pulses associated with the long-time behavior. For de- 
viations away from the quantization of the pulses in the 
case of a modulation of barrier profile, we found that, as 
expected, this led to the reappearance of the equilibrium 
noise contribution, which increases logarithmically with 
barrier opening time (see Fig. [T]). 

This work is supported by EPSRC-GB (Contract No. 
EP/D065135/1). 



APPENDIX A: NUMERICAL 
DIAGONALIZATION OF 



Riemann-Hilbert approach at finite temperature. The 
scattering matrix S{t) in takes the form 



Sit) 



), forte [0,t,] 
- sm I cos f ' 



and equals the identity otherwise. Following the no- 
tation in Ref. [10, we introduce the matrix R{X) = 

'\ 0^ 



S'e'^'^S'te^*^-^, with L = 
function In^ reads 







The characteristic 



In x(A) = Tr[n In R\ + Tr[ln(l - n + nR) - n In i?] (Bl) 



We describe the numerical procedure for diagonalizing 
n°"* in ([6|). We need a discretizcd finite dimensional ex- 
pression for XnX^ , where X stands for the transmission 
and reflection matrices A and B in frequency space, n 
is an infinite dimensional diagonal matrix with elements 
'^mn = '5m„^(— m + e), wlicrc e < o^q, to avoid to = in 
the 0-function. The matrix X is blockwise tridiagonal: 



X 



X2 Xo Xi 



V 



where Xo^i,2 is square matrix with truncated dimension 
M (see main text). Xi (X2) is lower (upper) triangular 
matrix with null diagonal elements. 

Exploiting the property B{t)B*{t) -f A(t)A*{t) 1, 
which in frequency space states BB^^ H- AA^^ = I (I is 
the infinite dimensional identity matrix), we find the fol- 
lowing properties for the finite block sub-matrices with 
dimension Af: X^LoC^iBI + A,A| 



I, B2BT 



AqAI = 0, and AjAJ, BiB 



BoBl 



0. 



Compute n°"* in ^ in frequency space and remove the 
part equivalent to n*" at low and high energies, we arrive 
at a much simplified form of n°"* suitable for numerical 
diagonalization 



BiBI + AiAI BiBJ + AiAJ 



BoBl-fAoAl I-B2B^-A2A^ 



The sought-after spectrum, rij, is obtained by diagonal- 
izing (|XT1) directly. 



where Tr operates on both time and channel space, n is 
the fermionic operator at finite temperature 1//?: 



n{t,t') 



V/3 



2tt sinh(7r(t-t')//3 + «0) 



Since there is no average charge transfer, the first term 
of dni]), Tr[7ilni?] = iX(Q), is equal to zero. The ECS is 
given by the second term. 



Inx(A) 



2n 



dX / dt tr 



dlnY: dluR 



dt dX 



where tr is a trace over channel space only. The Riemann- 
Hilbert solution Y{z) is a bounded matrix-valued ana- 
lytic function in the strip — /3/2 < Imz < f3/2 except 
along the cut z G [0,</], where Y_{t)Y^^{t) 
Y+{t) = Y{t+iQ) and Y_(t) = Y{t-iO) (see Ref. | 
Making the substitution sin ^ = sin ^ sin ^, R{X) can be 
diagonalized in a time- independent basis with eigenvalues 
g±iA, ^ Ugjjjg explicitly the finite temperature solution of 
the RH problem^ 



t • (Al) n(t) 



exp 



1 

27ri 



^ , , cosh(7rt//3) 7rlni?//3 



dt 



cosh(7ri7/3) siiih(7r(t - t')//3) 



we obtain 



APPENDIX B: RIEMANN-HILBERT SOLUTION 
TO THE FCS FOR BARRIER OPENING AT 
NON-ZERO TEMPERATURE 



Here we give a brief derivation of the FCS for the 
pure opening problem (region 2, Sec. Ill Cp within the 



l-x(A) = -^ln(^-j^^-p^j, 

with ^ as a cut-off. Taking the zero temperature limit 
(3 00, wc arrive at the characteristic function (j23p 
along with definitions (pi)) and . 
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